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ABSTRACT 



We introduce a new parameter A£ — the difference in magnitude between the 
red giant branch (RGB) bump and the point on the main sequence (MS) at the 
same color as the bump (which we call the "benchmark") — to estimate the he- 
lium content in old stellar systems. The A£ parameter is a helium indicator since 
an increase in helium makes, at fixed age and iron abundance, the RGB bump 
brighter and the MS benchmark fainter. Moreover, its sensitivity to helium is 
linear over the entire metallicity range. A£ is also minimally affected by changes 
of a few Gyr in cluster age, by uncertainties in the photometric zero-point, by 
the amount of reddening, or by the effects of evolution on the horizontal branch. 
The two main drawbacks of the A£ parameter include the need for precise and 
large photometric data sets from the RGB bump down to the MS benchmark, 
and a strong dependence of the AY/A£ slope on metallicity. To provide an 
empirical basis for the A£ parameter we selected almost two dozen relatively 
bright Galactic Globular Clusters (GGCs) with low foreground reddening, and a 
broad range of iron abundance (—2.45 <[Fe/H]< —0.70 dex). Moreover, the se- 
lected GGCs have precise, relatively deep, and homogeneous multi-band (BVI) 
photometry. We found that the observed A£ parameters and those predicted 
from a-enhanced evolutionary models agree reasonably well if we assume a pri- 
mordial helium content of Y=0.20 (abundance by mass). The discrepancy in 
the photometric B band becomes of the order of 4a (AB = 0.20 mag) only in 
the metal-poor regime. Comparison with evolutionary prescriptions based on a 
canonical primordial helium content (Y=0.245, AY/AZ [helium-to-metal enrich- 
ment ratio] =1.4) indicates that the observed A£ values are systematically smaller 
than predicted. This discrepancy ranges from 5a (AB = 0.26 mag) in the metal- 
rich regime to 10a (AB = 0.51 mag) in the metal-poor regime. The outcome 
is the same if predicted A£ parameters are based on evolutionary models with 
CNO enhancements in addition to a enhancements. The discrepancy becomes 
even larger if we extend the comparison to He-enhanced models. These findings 
support previous results (Meissner & Weiss 2006, A&A, 456, 1085; Di Cecco et al. 
2010, ApJ, 712, 527; Cassisi et al. 2011, A&A, 527, A59) suggesting that current 
stellar evolutionary models overestimate the luminosity of the RGB bump. We 
also found that including envelope overshooting can eliminate the discrepancy, 
as originally suggested by Alongi et al. (1993, A&AS, 97, 851); atomic diffusion 
and mass loss play smaller roles. The A£ parameter of GGCs, in spite of the 
possible limitations concerning the input physics of current evolutionary models, 
provides an independent detection of pre-stellar helium at least at the 5a level. 
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Subject headings: globular clusters: general - stars: evolution - stars: horizontal- 
branch - stars: red giant branch - stars: Population II 



1. Introduction 



The notion of possible variations in heliu m abundance between different Globular Clus- 



ters (GCs) dates back more than forty years (Ivan den Berghlll967f). in connect ion with the 
second parameter phenomenon ( VandenBerg 200ol ; Recio-Blanco et al. 2006 ). Recently, 



the possibility of different helium abundances among stars within individual GCs has also 
been suggested to explain not only the presence of multiple clumps and extended blue 
tails on t he horizontal branch, b ut also multiple parallel main sequences and subgiant 



branches (ID'Antona &: Caloil 120081 . and references therein). This working hypothesis relies 



on t he well establ i shed anti-co rrelations between the molecular band-strengths of CN and 



CH (ISmith 1 119871 ; iKraft I Il994l ) and between O-Na and Mg-Al measured in both evolved 



of GCs investigated with high-resolution spectroscopy ( 


Pilachowski, Sneden & Wallerstein 


1983; 


Ramirez & Cohen 


2002; 


Gratton, Sneden & Carretta 


2004). Furthermore, deep Hub- 



ble Space Telescope photometry has disclosed the presence of multiple stellar p opulations 
in several massive GCs. Together with the long-known case of u Centauri ( [Anderson 



GCs: NGC2808 (D'Antona et al. 



20021 ; iBedin et al. I l2004f) . parallel stellar sequences have now been detected in a numb er of 



2005 



NGC 1851 ( jCalamida et al. Il2007r iMilone et al. II2008I ). In some of these cases the multiple 



Piotto et al. 1 120071). M54 flSiegel et al. 1120071 ) and 



sequences might be explained by sub-populations dif f ering either in He a bundance ( INorris 



2004J : ID'Antona et al. 1 120051 : ID'Antona fc Caloil 120081 : iPiotto et al. II2007D . or in CNO abun 



dance. 

We still lack a direct detection of helium variations among the sub-populations in a GC, 
because helium absorption li nes in the vis ib le spectral region app ear only at effective temper- 
atures hotter than 10,000 K <jBehxll2003al fbL I Moehler et al. 1120041 ) . Such high temperatures in 
a GC are typical of hot and extreme HB stars but, unfortunately, their helium abundances 
cannot be adopted to constrain the original helium content because gravitational settling 
and /or radiative levitation can affect their surface abundances. In a recent investigation 
Villanova et al. I (120091 ) identified helium absorption lines in a few warm cluster HB stars 
that should be weakly affected by these mechanisms, but the range in effective temperature 
remains very narrow. 
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Interestingly enough, iPrestonl (120091 ) has identified both Hel and Hell emission and ab- 
sorption lines in field RR Lyrae stars. The lines appear during the rising branch of the light 
curve and appear to be the consequence of shock s propagating through the extended atmo - 



sphere soon after the phase of minimum radius (IBono et al 



994 



Chadid fc Gilletill996h . 



The idea is not new and dates back at least to IWallersteinl ( 119591 ). who detected helium 
emission lines in his seminal investigation of Type II Cepheids. The key advantage of RR 
Lyrae stars, when compared with similar detections of helium lines in hot and extreme HB 



stars (jBehrl l2003bl ) is that they have an extensive convective envelope, an d therefore they 



are m inimally affected by gravitational settling and/or radiative levitation ( IMichaud et al. 



20041 ). The drawback is that quantitative measurement of the helium abundance requires 
hydro dynamical atmosphere models simultaneously accounting for time-dependent convec- 
tive transport and radiative transfer together with the formation and the propagation of 
sonic shocks. 



A new approach has been suggested by iDupree et al. I (120111 ) , who detected the chro- 
mospheric Hel line at 10830 A in a few RG stars having similar magnitudes and colors in 
w Cen. The He line was detected in more metal-rich stars and it seems to be correlated with 
Al and Na, but no clear correlation was found with Fe abundance. However, a more detailed 
non-LTE analysis of the absolute abundance of He is required before firm conclusions can 
be drawn concerning the occurrence of a spread among the different sub-populations. 



To overcome the thorny problems connected with the spectroscopic measurements, llben 



( 119681 ) suggested the use of the R parameter — the ratio of the number of HB stars to the 
number of RGB stars brighter than the flat part of the HB — to estimate the initial He 
abundance of cluster star s. Two other indepen dent constraints on the overall helium content 
in GCs were suggested by lCaputo et al. I ( 119831 ). The A parameter, which is the difference in 
magnitude between the MS at (B-V) = 0.7 and the luminosity level of the flat part of the 
HB, and the A parameter, which is the mass-to-luminosity ratio of RR Lyrae stars . The pros 



and co ns of these parameters have been discussed in a thorough investigation by ISandquist 



( 120001 ). Theoretical an d empirical lim i tation s affecting the precisi o n of the R parameter have 
been also discussed by IZoccali et al. I ( 120001 ) , iRiello et al. I (120031 ) , ICassisi et al. I ( 120031 ) and 



Salaris et al. I (120041 ) 



The key feature of all three parameters is that they are directly or indirectly connected 
with the HB luminosity level. In spite of recent improvements in photometric precision and 
the increased use of a multi-wavelength approach, we still lack firm empirical methods to 
predict the HB luminosity level in GCs. This problem is partially due to substantial change 
in the HB morphology when moving from metal-poor (generally blue HB) to metal-rich (red 
HB) GCs. Additionally, we still lack a robust diagnostic to constrain the actual off-ZAHB 
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evolution of HB stars (IFerraro et al. Ill999t iDi Cecco et al. Il2010l ; ICassisi et al. Il201lh . 



To overcome these longstanding problems we propose a new parameter — which we call 
A£ — to estimate the helium content of GCs. The A£ parameter is the difference in magnitude 
between the RGB bump and a benchmark defined as the point on the (essentially unevolved) 
MS at the same color as the RGB bump. In the following we discuss the theoretical framework 
adopted and the sample of globular clusters we have selected to calibrate this new parameter. 



2. Theoretical framework 



To provide a quantitative theoretical framework fo r inferring the hel i um content of 
GCs we adopted the evolutionary models provided by Pietrinferni et al. ( 2004 . 2006 ) 1 ! 
These models were comp uted using a recent version of the FRANEC evolutionary code 
(jChieffi fc Stranierd Il989l ). The reader interested in a detailed discussion of the adopted 
input physics is referred to the above papers. Her e we mention only that we adopted the 



radiative opacity tables^ from the Live rmore group 



peratures higher than 10 4 K, and from lFerguson et al- 



exias &: Rogerslll996f ) for interior tem- 
(2005) for the atm ospheres. Thermal 



conduction is accounted for following the prescriptions by lPotekhinl (ll999l). We have update d 
the plasma-neutrino energy loss rates according to the results provided bv lHaft et al. (1994). 



The nuclear reaction rates have been updated using the NACRE data base ( lAngulo et al. 



19991 ). with the excepti on of the 12 (?(a;7) 16 reaction; for this reaction, we adopted the 
recent determination by IfCunz et al. ( 2002 ). We also adopted the equation of state (EOS) 
provided by A. Irwiro 



The adopted a-enha nced mixture follows the prescriptions bv lSalaris et al. I (119971) an d 
by Salaris fc Weiss ( 1998 ). with a Solar ir on content of 0.0198 from Grevesse &: Noelsl ( 1993 ). 
a primordial helium c ontent of Y=0.245 (jCassisi et al. 1120031) . and a helium-to- metal enrich- 
ment of AY/AZ=1.4 ( IPietrinferni et al. 1120041 ) . The cluster isochrones were transformed into 
the observational plane using the bolometric co r rectio ns (BCs) and the color-temperature 
relations (CTRs) provided by lCastelli &: Kuruczl (120031 ). 



Note that together with the a-enhanced models, we al so considered cluster iso chrones 
constructed assuming CNO-Na abundance anti-correlations (IPietrinferni et al. II2009I ). These 



1 See the BaSTI web site at the following URL: |http : //193 . 204 . 1 . 62/indexT html 

2 The reader interested in a detailed description is referred to the the following URL: 



http: / /opalop acity.llnl.gov/opal.html| 

3 A detailed description of this EOS can be found at the following URL: http://freeeos.sourceforge.net 
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sets of isochrones have the advantage of having been constructed for the same iron abun- 
dances as the a-enhanced models. The heavy- element mixture adopted in these evolutionary 
models is defined as extreme (see Table 1 in iPietrinferni et al. 1 120091 ). and when compared 
with the canonical a-enhanced mixture includes an increase in N and Na abundance of 1.8 
and 0.8 dex and a decrease in C and O abundance of 0.6 and 0.8 dex, respectively. We did 
not include t he anti-correlation b etween Mg and Al, since its impact on evolutionary models 
is negligible ( ISalaris et al. 1 120061 ) . The transformation of these evolutionary models into the 
observational plane does require a set of bolometric corrections (BCs) and color-temperature 
relations (CTRs) comp uted for the same mixtures; these be came available only very recently 
( ISbordone et al. 1120101 ) . However, IPietrinferni et al. I (120091 ) found that BCs and CTRs com- 



puted for simple a-enhanced mixtures mimic the same behavior. Moreover, |Pi Cecco et al. 



( 120101 ) found that BCs and CTRs hardly depend — at fixed total metallicity — on changes in 
the detailed mixture. The total metallicity of the atmosphere models adopt e d to t ransform 
the a- and CNO-enhanced models was estimated using the ISalaris et al. I ( 119931 ) relation 
with [a/Fe]=0.62 dex. 

Finally, we mention that current isochrones together with similar predictions available 
in the literature (jVandenBerg et al. Il2000l ; iDotter et al. 1120071 ; iBertelli et al. 1120081 ) give ages 
for G Cs that are, within the errors, consistent w ith the WMAP age for the Universe of 13.7 



Gyr (IKomatsu et al 



2011 



the Galactic GCs (iMarin-Franch et al. 112009 : ICassisi et al. 1120111 ) 



Larson et al. 1120111). and provide a qu ite similar age ranking for 



3. Optical data sets 

To provide a solid empirical basis for validating the A£ parameter we selected a sam- 
ple of relatively luminous GCs (for good statistics), covering a broad range in metallicity 
(-2.45<[Fe/H]<-0.70), with different structural parameters, and all having relatively low 
foreground reddening (E(B-V) < 0.14, except for NGC6352 with E(B-V) = 0.21 mag). 
Moreover and even more importantly, we selected GCs for which precise and relatively deep 
photometry is available in the BVI photometric bands. To reduce the possibility of any 
subtle errors caused by the absolute photometric calibration or by the reduction strategy, we 
only adopted catalogs available in the data base maintained by co-author PBS. We ended 
up with a sample of almost two dozen GCs; their names, reddenings, metallicities and total 
visual magnitudes are listed in the first four columns of Ta ble 1. Note tha t the c luster iron 



abundances are based on the metallicity scale provided by Kraft fc Ivans! (120031 ) based on 



Fell lines in high-resolution spectra. For the clusters for which the Kraft & I vans metallicity 



was not available we adopted the iron abundances on the metallicity scale by lCarretta et al. 
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( 120091 ) calibrated using Fel lines in high-resolution spectra. The mean difference between the 
two metallicity scales is minimal (<0.1 dex) and vanishingly small in the metal-poor regime 
(see their Fig. 10). 

All o f these c l usters have been calibrated t o the Johnson- Cousi ns BVI photometric 
system of lLandoltl (119731 Il992l ) as described by IStetsonl (120001 . 120051 ). In every case the 
corpus of observations is quite rich, ranging from a minimum of 118 (NGC6362 in B) to a 
maximum of 3939 (NGC 6341 in /) individual CCD images per filter in each cluster (median: 
946 images per filter per cluster). The instrumental observations have been referred to the 
standard system by means of local standards within each field. The GCs in our sample have 
calibrated data in all three of the B, V, and / bands based on a minimum of 118 (NGC 6362) 
and a maximum of 3690 (NGC 6341) local standards (median: 972). Observations for each of 
the different clusters were obtained during a minimum of 7 (NGC 6541) and a maximum of 44 
(NGC 6341) different observing runs (median: 16). It must be stressed, however, that since 
our proposed methodology relies on the difference in apparent magnitude between two stellar 
sequences at fixed color absolute photometric zero points are unnecessary, and possible errors 
in the color transformations of the photometric calibration are irrelevant (provided any such 
errors are unaffected by stellar surface gravity). The principal observational requirement 
of this work is that the photometry be consistently precise and linear over the required 
magnitude range. We believe that our data set, based on a very large number of individual 
CCD images from many different observing runs and a substantial and homogeneous body 
of local photometric standards for each cluster, guarantees that this requirement is met to a 
higher degree than in any previous photometric surveys of Milky Way GCs. 

For each GC we estimated the ridge line in both the V,B-V and the V, V-I color- 
magnitude diagram (CMD). The exact approach adopted to estimate the ridge lines will 
be described in a future paper (Ferraro et al. 2011, in preparation). Fig. 1 shows three 
different CMDs for three clusters that are representative of metal-rich (47 Tuc), metal- 
intermediate (M3) and metal-poor (M92) GCs. The red solid curves display the ridge lines 
that we have estimated. The approach we adopted to estimate the magnit ude and the color 



of the RGB bump has already been described in |Pi Cecco et al. I (120101). The differen ce 



between the current \^-band RGB bumps and those estimated by iDi Cecco et al. I (J2010f) is 
on average 0.00±0.04 m ag (12 GCs). The difference compared to the RGB bumps estimated 
bv iRiello et~al~l J2003I) is +0.01±0.07 mag (12 GCs), while the difference with the RGB 
bumps provided by lCassisi et al. I ( 120111 ) is -0.03±0.05 mag. We performed a series of tests to 
estimate the accuracy of the position of the bump along the RGB using different colors (B-I, 
B-V, V-I) and we found that among them V-I provides the most precise measurements. We 
plan to provide a more detailed discussion concerning the pros and cons of both optical and 
near-infrared colors to estimate the A£-parameter in a future paper. The B, V magnitudes 
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012301230123 
B-I [mag] B-I [mag] B-I [mag] 



Fig. 1. — Top - B,B-I color-magnitude diagram from left to right of a metal-rich ([Fe/H] =- 
0.70, 47 Tuc), a metal-intermediate ([Fe/H] =-1.50, M3) and a metal-poor ([Fe/H] =-2.38, 
M92) GC. The red solid curves show the cluster ridge lines. The red arrows mark the position 
of the RGB bump, while the green asterisks show the MS benchmark. The vertical red lines 
show the difference in magnitude between the RGB bump and the MS benchmark. Middle 
- Same as the top, but in the V,B-I color-magnitude diagram. Bottom - Same as the top, 
but in the I, B-I color-magnitude diagram. 
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and the V-I colors of the RGB bumps are listed in columns 5, 6 and 7 of Table I. The 
B, V magnitudes of the MS benchmark having the same V-I color as the RGB bump were 
estimated by interpolating the ridge line (see the green asterisks plotted in Fig. 1) and they 
are listed in columns 8 and 9 of Table 1. 

The data in Table 1 confirm the visual impression from Fig. 1 that the luminosity of 
the RGB bump varies inversely with the metallicity of the cluster. The A£ parameter is 
therefore more prone to Poisson statistics and random photometric errors at the metal-poor 
end since in this regime the bump falls in a region of the RGB where the lifetime is shorter 
than at fainter magnitudes and the relative number of RG stars steadily decreases. The 
color of the bump also becomes redder relative to the turnoff due to the slope of the RGB, 
and the apparent magnitude of the MS benchmark becomes correspondingly fainter. 

4. Comparison between theory and observations 

We estimated the A£ parameter in two different bands (B and V); these results are listed 
in columns 10 and 11 of Table 1, together with their intrinsic errors, and plotted in Fig. 2. 
The adopted cluster isochrones were retrieved from the BASTI data base. In particular, 
we estimated the predicted A£ parameter using cluster isochrones based on evolutionary 
models constructed for iron abundances covering the entire range typical of old, low-mass 
stars, a-enhanced ([a/Fe]=+0.40 dex) mixtures and a mass-loss rate d la Reimers with 
77=0.4. We adopted five different helium abundances, namely Y=0.20, 0.245, 0.30, 0.35 and 
0.40 (relative abundance in terms of mass fraction). The models with Y=0.20, 0.30, 0.35 and 
0.40 were constructed assuming a constant helium content over the entire metallicity range, 
while those with Y=0.245 were constructed assuming a helium-to-metal enrichment ratio 
AY/AZ =1.4 (see §2). Note that the models with Y=0.20 have been specifically computed 
for this experiment, but they will become available at the web site mentioned above. To 
constrain the possible dependence on cluster age, we adopted three different values spanning 
the generally accepted range of GC ages, namely 10, 12, 14 Gyr (see dashed, dotted and 
solid lines in Fig. 2). 

The predictions for Y=0.40 were not included in Fig. 2 because the bump was identifiable 
only in a few metal-rich compositions and for the older ages. The increase in helium causes, 
at fixed age and metal content, a decrease in the envelope opacity. This means that the 
maximum depth reached by the outer convective region during the first dredge-up becomes 
shallower and the discontinuity in the hydrogen profile left over by the retreating convective 
envelope is milder. Therefore, the H-burning shell crosses the chemical discontinuity at 
brighter magnitudes during faster evolutionary phases and, as a result, the bump becomes 
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Fig. 2. — Top - Difference in the 5-band magnitude between the RGB bump and the 
MS benchmark for the selected GCs as a function of the iron abundance. The selected 
GCs are marked with different symbols and colors. The dashed, dotted and solid lines 
display theoretical predictions for three different cluster ages (10,12,14 Gyr) according to 
the BASTI data base. The adopted cluster isochrones are based on evolutionary models 
computed assuming an a-enhanced mixture, different iron abundances (-2.6< [Fe/H] < -0.3 
dex) and five different He contents (Y=0.20-0.40). Predictions for Y=0.40 are not plotted, 
since the bump was identifiable only in a few isochrones. Bottom - Same as the top, but for 
the l/-band magnitude. 
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less and less evident. 

The main advantages in using the A£ parameter to estimate the helium content are 
the following: i) An increase in He content causes, at fixed cluster age and iron content, a 
decrease in the envelope opacity and an increase in the mean molecular weight. This means 
that an increase in He makes the RGB bump systematically brighter and the MS benchmark 
systematically fainter. The sensitivity to helium is very good, and indeed the derivatives 
for a-enhanced structures, at fixed iron and age, attain values similar to the A parameter 
(see Table 2). it) The sensitivity to age is small, and indeed the derivatives for a-enhanced 
structures, at fixed iron and helium, attain minimal values. Hi) There is no dependence on 
photometric zero-point or interstellar extinction or reddening, iv) The derivative is almost 
linear over the entire metallicity range (see Fig. 2). This evidence is supported by both theory 
and observation. The R parameter, in co ntrast, shows a well defined change in the slope 



for [Fe/H]fti-0.70 dex (IZoccali et al. 1 120001 ) . v) There is no dependence on HB evolutionary 
effects. In contrast, the R, A, and A parameters are all affected by the off-ZAHB evolution 
of HB stars. 

The use of the A£ parameter has two main drawbacks: i) Precise, linear photometry 
from the RGB bump down to the MS benchmark is required for photometric data sets large 
enough to properly identify the RGB bump, ii) There is a definite metallicity dependence, 
and indeed the derivative (~1 mag/dex) is at least 0.7 dex larger than for the other helium 
indicators (see Table 2). The above points confirm that robust He estimates based on the 
A£ parameter do require precise photometry and accurate spectroscopic iron abundances. 

The positive features of the predicted A£ parameter allow us to derive an analytical 
relation with iron and helium content. Using a multi-linear regression over the entire set of 
A£ B - and A£ v parameter, we found: 

Y = -0.525 ± 0.037 + 0.158 ± 0.007 A£ B + 0.190 ± 0.008 [Fe/H] a = 0.018 

[1] 

Y = -0.537 ± 0.037 + 0.157 ± 0.007 A£ y + 0.180 ± 0.008 [Fe/H] o = 0.018 



where a is the standard deviation and the other symbols have their usual meaning. Note 
that we have neglected the age dependence because the dependence on this parameter is 
imperceptible within the typical age range of GCs. 

The above analytical relations do allow us to provide individual helium estimates for 
GCs. However, the comparison between theory and observation indicates that the pre- 
dicted A£ parameters are systematically larger than observed and the discrepancy steadily 
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increases when moving from lower to higher helium contents. Moreover, the discrepancy, at 
fixed helium content, increases when moving from metal-rich to metal-poor GCs. To repre- 
sent the observed, empirical values of A£ as a function of metallicity, we performed linear 
least-squares fits to the data in Table 1 and obtained: 

A£ B = 4.75 ± 0.13 - 1.04 ± 0.08 [Fe/H] a = 0.08 mag 

[2] 

A£ v = 4.85 ± 0.12 - 1.00 ± 0.07 [Fe/H] a = 0.09 mag 

The above fits were performed by using as weights the intrinsic errors on the A£ values. 
The anonymous referee noted that the slope of the A£ relations might not be linear over the 
entire metallicity range. In particular, the slope seems to become steeper in the metal-rich 
regime. However, this inference is based on only two GCs (NGC 6352, NGC104). Additional 
data in the high metallicity tail of GGCs are required to reach more quantitative conclusions. 
However, to constrain the impact of a possible change in the slope we performed the same 
fits by using only GCs more metal-poor than [Fe/H] <-l dex. We ended up with a sample 
of 20 GCs and found: 

A£ B = 5.00 ± 0.15 - 0.90 ± 0.09 [Fe/H] a = 0.13 mag 

[3] 

A£ v = 5.18 ± 0.13 - 0.81 ± 0.07 [Fe/H] a = 0.17 mag 

The new fits are characterized by slopes that are slightly shallower and zero-points that 
are mildly larger, but the two sets of least-squares solutions agree within la. 

The small intrinsic dispersions of the above empirical relations suggest that the A£ 
parameter can be adopted to identify GCs with peculiar He abundances. The A£ values 
listed in Table 1 and the analytical relations based on the less metal-rich GCs indicate that 
the possible dispersion in helium among the target GCs is minimal. We found that only 
two GCs — M13 and NGC 6541 — show differences larger than 1.5a both in the A£ B and 
A£ v parameters. However, the discrepancy concerning NGC 6541 should be treated with 
caution, since it and NGC 6352 are the two GCs in our sample with the largest reddening. 
We are only left with M13 for which the discrepancy ranges from 1.3 to 2.6<r in the V and 
in the B bands, respectively. By using the theoretical relations [1] we found that the quoted 
differences imply on average a difference in helium content, when compared with the bulk of 
GCs, that is of the order of la (AF=0.031±0.025). The lack of a clear evidence of helium 



13 



enhancement in M13 soundly supports the recent results obtained by lSandquist et al.l (120001 ) 
using the R parameter. Unfortunately, we cannot extend the current differential analysis 
concerning the helium content to more metal-rich GCs, due to the scarcity in our sample of 
GCs in this metallicity regime. 

In this context it is also noteworthy that if the He content is the second parameter 
affecting the HB morphology (ISuda fc Fujimotol 120061 ). the A£ parameter should be able 
to prove it. At the present moment, the empirical evidence does not allow us to reach a 



firm conclusion. The classical couple of second parameter GCs M3 and M 13 (I Johnson et al. 
20051 ) — i.e. two GCs that within the errors have the same chemical composition, but different 
HB morphologies— show a mild difference («2<r) in A£ (A£ B = 6.38±0.21 [M3] vs 6.74±0.09 
[M13]; A£ v = 6.43±0.07 [M3] vs 6.73±0.15 [M13]) level. On the other hand, the other 



classical couple of second parameter GCs NGC 288 and NGC 362 (jBellazzini et al. 1120011 ) 
do not show, within the errors, any difference. The current scenario might be confused 
by systematic uncertainties affecting either the chemical composition (iron, a-enhancement, 
CNO-enhancement) or the cluster age which are not directly correlated with a difference 
in He. Clearly, the simple-minded notion that He is the second parameter is not obviously 
correct; other factors, which are not immediately obvious, will have to be taken into account. 

The comparison between the empirical and predicted relations indicates that the dis- 
crepancy in the B band for the lowest helium content (Y=0.20) is negligible in the metal-rich 
regime ([Fe/H] ~-0.8 dex) and slightly smaller than 5a (AB = 0.26 mag) in the metal-poor 
regime ([Fe/H] dex). The difference in magnitude with canonical helium models ranges 

from AB = 0.26 (5a, [Fe/H] «-0.8 dex) to AB = 0.51 (10a, [Fe/H] w-2.2 dex) and becomes 
even larger for He-enhanced models (Y=0.30: 0.64, 0.89 mag; Y=0.35: 0.95, 1.21 mag in the 
metal-rich and in the metal-poor regime, respectively). The difference hardly depends on the 
adopted photometric band, and indeed the V-band difference with the lowest helium models 
is negligible in the metal-rich regime and equal to AV^ = 0.20 in the metal-poor regime (4a, 
[Fe/H] ps-2.2 dex). The difference with canonical helium models ranges from AV = 0.25 
(ha, [Fe/H] w-0.8 dex) to AV = 0.45 (9a, [Fe/H] «-2.2 dex). 

The above results indicate either that the primordial h elium content is roughly 20% lower 
than estimated using t he CMB (WMAP Isteigmanlboioh and the HI I regions in metal-poor 
blue compact galaxies (jOlive fc Skillmanl 12004 ; llzotov fc Thuanll2010l ; iPeimbert et al. 1120 lOh 
or that theory overestimates the luminosity of the RGB bump . The latter working hypoth- 
esis is also supported by recen t results by iMeissner fc Weissl ( 120061 ) and more recently by 
Di Cecco et al. I (120101 ) and by ICassisi et al. I (120111 ) . The luminosity predicted by current 
evolutionary models is brighter than observed in actual GCs. Note that the current dis- 
crepancy is more striking than the previous estimates, since the A£ parameter is minimally 
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Fig. 3. — Same as Fig. 2, but the stellar isochrones are based on evolutionary models 
constructed assuming an a- and CNO-enhanced mixture. The adopted He abundances range 
from Y=0.24 to Y=0.40. 
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affected by uncertainties in distance, age, reddening, and HB luminosity level. 

To further constrain the nature of the possible culprits, we also tested the role played 
by possible overabundances in both a- and CNO elements simultaneously. The comparison 
of the observations with predictions based on evolutionary models constructed assuming 
an extreme chemical mixture (see §2) is shown in Fig. 3. Note that the comparison was 
performed at fixed iron content, since the iron abundance is found to be constant even in 
GCs showing well defined anti-correlations in the lighter metals. The data plotted in this 
figure indicate that the impact of a- and CNO-enhanced mixture on the current discrepancy 
is small. 



4.1. Possible stellar modeling venues 



Some crucial details of both micro- and macro-physics assumptions in the modeling of 
low- mass stars are far from bein g settled. Som e of th em may affect, directly or indirectly, 
the bump position in isochrones. ICassisi et al. I (120111 ) already discussed a number of them 
(their Sect. 4.1), either estimating their impact on the bump and on the turnoff position 
from their own calculations or using earlier results available in the literature. Among these 
uncertainties atomic diffusion and convective overshooting are the most interesting ones. 
Here we give an example of their possible influence on evolutionary tracks in the theoretical 
H-R diagram for a model with M = 0.8 M & and chemical composition Z=0.001, Y=0.250. 
The adopted chemical composition is a-element enhanced: it corresponds to [Fe/H] = -1.56 
dex. Such a star has a turn-off age of ~ 13 Gyr, and therefore its RGB is quite representative 
for the RGB of an is ochrone of this age. T he evolutionary tracks were constructed using 
the GARSTEC code ( Weiss fc Schlattl 2008T). and the evo lutionary models we re computed 



with the updated 14 N(p, 7) 15 reaction (iMarta et al. 1120081 ). which was shown (I Weiss et al. 
20051 ) to increase the bump luminosity by ~ 0.06 mag. Note that the BASTI evolutionary 



tracks do not account for this reaction; its inclusion would increase the discrepancy between 
the predicted and the observed A£ parameter by ~0.15 mag. 

We start with a canonical model without convective overshooting and neglecting atomic 
diffusion (black line, top panel of Fig. 4)). In this case, the bump (which is not very pro- 
nounced for these stellar parameters) is located almost exactly at log L/L Q = 2.0. The 
inclusion of atomic diffusion (red line) decreases — as it is well-known — the main-sequence 
lifetime by ~0.5 Gyr, but has only a minimal impact on the bump luminosity (triangle). In- 
deed, it becomes slightly brighter by only Ri0.03 mag (~0.01 d ex in log L). We then in cluded 
the overshooting, treated as a diffusive process as described in lWeiss fc Schlattl! (120081 ). with 
a scale-length parameter of 0.04. The adopted value is approximately twice the value needed 
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Fig. 4. — Top - Zoom of the evolutionary tracks around the position of the RGB bump for a 
model with M = 0.8 M Q and chemical composition Z=0.001, Y=0.250, constructed assuming 
different assumptions concerning atomic diffusion, envelope overshooting and mass loss. The 
canonical evolutionary track (black line) neglects the above mechanisms. The position of 
the RGB bump is marked with different symbols. See text for more details. Bottom - Same 
as the top, but as a function of the logarithmic effective temperature. 
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to fit the CMDs of open clusters (see, e.g. iMagic et al. 1 120 10h . Convective overshooting 



is applied at all the convective boundaries, and in particular at the lower boundary of the 
convective envelope. However, it is suppressed during the transient convective core phase 
of the pre-main sequence evolution. The blue line plotted in the top panel of Fig. 4 shows 
that inclusion of the envelope overshooting makes the bump (filled blue square) fainter by 
0.38 mag (0.15 dex in logL) when compared with the canonical evolutionary track. This 
is more than enough to remove the discrepancy between the predicted and observed A£ 
parameter. Note that the difference in luminosity between the bump in the canonical evolu- 
tionary track and in the track including overshooting is causing a minimal change in effective 
temperature (A T e ff ~65 K, see the bottom panel of Fig. 4), due to the steep slope of the 
RG B along these evolu tionary phases. This finding confirms the results originally suggested 



by lAlongi et al. I (Il993[ ) concerning the dependence of the lum inosity of the RGB bump on 



the envelope overshooting (see also ICassisi et al. 1 12002k l201lh 



As a further test, we included atomic diffusion together with envelope overshooting. 
The green line plotted in Fig. 4 shows that the bump (open green square) is getting even 
fainter, but the difference compared to the track including only the envelope overshooting 
is ~0.05 mag (0.02 dex in logL). Finally, together with envelope overshooting and atomic 
diffusion we also included mass loss according to the Reimers formula with the free parameter 
1] = 0.6 (cyan line). Note that this value is larger than the value typically adopted, rj = 0.4. 
However, the impact on the bump position (cyan filled circle) is negligible, the difference 
being of the order of a few hundredths of a magnitude. 

The inclusion of diffusion and / or convective overshooting has a minimal influence on the 
MS benchmark, since at the effective temperature of the bump, the corresponding mass on 
the main sequence is hardly affected (see §4.1). Therefore, the luminosity of t he MS bench- 



mark at fixed T e fj is unaffected. This is the main difference with the results by lCassisi et al. 



( 120111 ). since along the MS they adopted the turn-off as reference point. Our indicator is 
age-independent to a high degree of confidence. However, there is an indirect effect: the 
convective overshooting moves the bump downwards along the RGB, which implies a blue- 
ward change of the bump logT e g, in this case by +0.0055 dex. The MS benchmark has 
therefore to move hotter (and brighter) on the MS by the same amount. Given the MS 
slope of A log Lj A log T c g = 7.5, this implies a luminosity increase of ~0.1 mag (~10% in 
luminosity), thus decreasing the discrepancy between theory and observation even more. In 
summary, our test calculations show that reasonable assumptions concerning the convec- 
tive overshooting from the bottom of the convective envelope goes in the right direction to 
significantly reduce the discrepancy, with the atomic diffusion playing a secondary role. 

It is worth noting that the quoted numerical experiments were performed at fixed chem- 



-18 - 



ical composition and cluster age. This means that the possible explanation we are suggesting 
to remove the discrepancy between predicted and observed A£ parameters needs to be tested, 
together with the new 14 N(p, 7) 15 reaction, in both the metal-poor and the metal-rich 
regime (Weiss et al. in preparation). 



Finally, it is worth mentioning that I Charbonnel fc Zahnl ( 120071 ) , to explain the abun- 
dance anomalies of RGB stars, suggested the occurrence of a thermohaline instability along 
the RGB when the hydrogen burning shell reaches the layers with uniform chemical compo- 
sition left behind by the first dredge-up. However, the occurrence of this phenomenon has 
a minimal impact on the luminosity of the RGB bump, since the hydrogen burning shell 
reaches the chemical discontinuity at earlier evolutionary phases. Moreover, recent evolu- 
tionary calculations, including state-of-the-art composition transport predi ctions, indicate 



that the thermohaline region does not reach the outer convective envelope (I Wachlin et al. 



201ll ). This limits the possibility for a non-canonical mixing process and for relevant changes 



in the chemical stratification above the hydrogen burning shell. 



5. Conclusions and final remarks 

We introduce a new A£ parameter — the difference in magnitude between the RGB bump 
and the MS benchmark (defined as the point on the MS with the same color as the RGB 
bump) — to estimate the helium content of globular clusters. The advantages of the A£ pa- 
rameter are the following: i) An increase in He content causes, at fixed age and iron content, 
an RGB bump systematically brighter and an MS benchmark systematically fainter. The 
sensitivity to He is very good, and indeed the derivatives for a-enhanced structures, at fixed 
iron and age, attain values (6.47-6.56 mag) similar to the A parameter (6.63 mag), ii) It has 
a minimal sensitivity to age, and indeed the derivatives for a-enhanced structures, at fixed 
iron and helium, are minimal (-0.03 mag/Gyr). in) It has no dependence on uncertainties 
in the photometric zero-point or the amount of reddening and extinction, iv) Its sensitivity 
to helium is almost linear at every metal abundance, vi) It is minimally affected by evolu- 
tionary effects. In contrast, the R, A, and A parameters are all affected by the post-ZAHB 
evolution of HB stars. 

The A£ parameter has two main drawbacks: i) Precise linear photometry from the RGB 
bump down to the MS benchmark is required, with large enough photometric data sets to 
properly characterize the RGB bump, ii) The slope of the A£ versus Y relation has a strong 
metallicity dependence. The derivative is ~-l mag/dex/dex, i.e., a factor of 0.7 larger than 
the A parameter and ~1 dex larger than the R and A parameters. 
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In passing we note that the A£ parameter, at constant helium-to-metal enrichment ratio, 
can also be adopted as a metallicity indicator. The slope of the metallicity dependence is 
similar to the slope of the RGB bump (1 vs 0.9), but the A£ parameter is minimally affected 
by uncertainties in cluster distance and in reddening. 

We have measured the A£ parameter in almost two dozen GGCs. The sample was 
selected according to the following criteria: homogeneous photometry; low reddening; a 
broad metallicity range (-2.5 < [Fe/H] ^-0.7 dex); and a range of structural properties. We 
found that, at first glance, theory and observation agree quite well if we assume a primordial 
He content of Y=0.20, and indeed they agree reasonably well in both the metal-rich and 
the metal-intermediate regime. However, the observed values in the metal-poor regime are 
5a (AB=0.26 mag) smaller than predicted. The difference with models assuming a more 
canonical primordial He abundance is larger, and ranges from 5a (AB = 0.26 mag) in the 
metal-rich regime to Aa (AB = 0.51 mag) in the metal-poor regime. The difference becomes 
even larger when moving to higher He contents (such as are assumed in some current inter- 
pretations of the multiple populations found in some GCs). These results support previous 
evidence tha t current evolutionary prescriptions generally overestimate the lumi nosity of the 
RGB bump rtMeissner & Weissll2006[ bi Cecco etaTlboiot ICassisi et alTlboilh . 



We computed a series of alternative evolutionary tracks at fixed age (t = 13 Gyr, 
M=0.80 M ) and chemical composition (Z=0.001, Y=0.250) and we found that plausible 
assumptions concerning the treatment of the envelope overshooting can account for the 
discrepancy between theory and observations. The atomic diffusion and the mass loss may 
also play secondary roles. For a tho rough discussion of th e possible impact of other physical 
parameters the reader is referred to ICassisi et al. I ( 120111 ) . 



The current debate concerning the primordial helium content is far from being settled. 
The primordial helium abundance mainly affects the third and the fourth acoustic peak 
of the CMB angular power spectrum. The seven years of observations by WMAP have 
provided an opportunity to improve the accuracy of the third peak. This means that He 
can be treated, for the first time, as a fitted paramete r in th e ACDM model. In particular, 
using a flat prior with 0.01 < Y P <0.8, lLarson et al~l (hoilh found Yp=0.28i[^5 (see their 
Table 8). This confirms the existe nce of pre-stella r heliu m at least at a 2a level. These 
findings were strongly supported by lKomatsu et al. I ( 120111 ) who com bined the seven years of 



WMA P data with data from small-scale CMB experiments (ACBAR, iReichardt et al. 112009 



QUaD. lBrown et al. II2009I ). They found a 3a detection of primordial He (Yp=0.326±0.075 
[68% CL], see their figures 10 and 11). To further constrain the primordial helium they 
adopted a very conservative upper limit (Yp<0.3) and they found (0.23 <Yp<0.3 68% CL). 
These geometric determinations of the primordial helium abundance do not by themselves 
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rule out an initial Y value as small as 0.20. 



According to the most recent calculations of state-of-the-art Big Bang Nucleosynthesis 
(SBBN) models the primordial helium content i s Yp 0.2483 ± 0.0005. The interested 



reader is referred to the comprehensive review by ISteigmanl ( 120071 ). This value agrees quite 
well with helium abundance estimates based on metal-poor extragalactic HII regions, namely 
Yp » 0.240 ± 0.006, but the impact on th e se measurements of possible systematic errors 
is sti l l controversial (Izotov fc Thuanl 12004 ; lOlive fc Skillmanl l2004t iFukugita fc Kawasaki 



2006; Peimbert et al. 



20071 ). 



In light of the above findings concerning the primordial He content, we followed this 
path: we estimated the primordial helium content using the analytical relations given in §3 
and the A£ parameters listed in Table 1. We found Y P = 0.181 ±0.014 (<r=0.034, A£ B ) and 
0.183±0.014 (a=0.033, A£ v ). The current Y P values were estimated as weighted means, the 
errors are the standard error of the mean and account for the uncertainties in the coefficients 
of the analytical relations, while the cr's are the standard deviations. It is worth noting 
that the Yp values we found are upper limits to the primordial helium content, since some 
extra helium is dredged up to the surface by the deep penetration of the convective envelope 
during the so-called first dredge-up phase. However, the amount of extra helium is very 
limited, and indeed at fixed cluster age (12 Gyr) it ranges from 0.008 ([Fe/H] =-2.62 dex, 
M[TO]=0.80 M ) to 0.02 ([Fe/H] =-0.70 dex, M[TO]=0.90 M Q ). 

The above estimates are within la of the esti mates provided by the C MB experiments, 
and by the analysis of lower main sequence stars (ICasagrande et al. 1120071 ). but are 2a from 
the estimates provided by the SBBN models and by the measurements of extragalactic HII 
regions. 

We also assumed that the current evolutionary models, when moving from the metal- 
poor to the metal-rich regime, on average underestimate the luminosity of the A£ parameter 
by 0.35 mag because of incomplete input physics. Therefore, we can apply — as a preliminary 
correction — a systematic shift of 0.35 mag to the analytical relations. Note that the coeffi- 
cients of the A£ parameter are of the order of 0.016, so the current shift causes an increase 
in helium of a few percent. Under these assumptions we found a primordial helium content 
of Y P =0.236±0.014 (A£ B ) and 0.238±0.014 (A£ y ). 

Consideration of the above results indicates that we still lack a convincing determination 
of the primordial He content based on stellar observables. However, if we trust the input 
physics currently adopted in evolutionary models, the existence of pre-stellar helium is estab- 
lished at least at the 5a level and increases to at least a 6a level, if we assume that current 
models underestimate the luminosity of the RGB bump by a few tenths of a magnitude. 
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No doubt the new and accurate CMB measurements by PLANCK (llchikawa et al. 



20081 ) will shed new light on the long-standing problem of the abundance of the second most 



abundant element in the Universe. The stellar path seems also very promising, but we still 
need to pave it and GAIA is a fundamental step in this direction. 
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Table 1. Apparent B- and K-band magnitude of the RGB bump and MS benchmark, together with the (V — I) 
color of the bump. The reddening, the mean metallicity and the total visual magnitude for each cluster are also listed. 
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a Cluster reddening according to lHarris 1996 ] l|2010 edition!) 

b Cluster absolute visual magnitude according to lHarris 1996 ] {2010 edition!) 

c Cluster iron abundance according to lKraft fc Ivan ] l|200j l2004h 

d Cluster iron abundance according to lCarretta et a~ ] j2009h 
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Table 2. Derivative as a function of iron, helium and cluster age for the most popular 
parameters adopted to estimate the helium content in GCs. 



Derivative 
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a The derivative (mag) was 


computed at fixed 


iron content and cluster a{ 


;es over 


the entire 


range 



of He abundances covered by a- and CNO-cnhanced models (BASTI data base). 

b The derivative (mag/dex/dex) was computed at fixed helium content and cluster ages over the 
entire range of iron abundances covered by a- and CNO-cnhanced models (BASTI data base). 

c The derivative (mag/Gyr) was computed at fixed helium and iron abundances over isochroncs 
with ages ranging from 10 to 14 Gyr (BASTIatabase). 

d The derivatives were computed adopting a-enhanced models (BASTI data base). 



